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, We demonstrate that the binding energies and widths of ?7-mesic nuclei 

0\ , 

. depend strongly on subthreshold ry-nucleon interaction. This strong depen- 

(N 

dence is made evident from comparing three different r/-nucleus optical poten- 
tials: (1) a microscopic optical potential taking into account the full effects of 
^ , off-shell r]N interaction; (2) a factorization approximation to the microscopic 

optical potential where a downward energy shift parameter is introduced to 
approximate the subthreshold r]N interaction; and (3) an optical potential 
using on-shell r]N scattering length as the interaction input. Our analysis 
indicates that the in-medium r]N interaction for bound-state formation is 
about 30 MeV below the free-space r/N threshold, which causes a substantial 
reduction of the attractive force between the r/ and nucleon with respect to 
that implied by the scattering length. Consequently, the scattering-length 
approach overpredicts the and caution must be exercised when these latter 
predictions are used as guide in searching for r/-nucleus bound states. We 
also show that final-state-interaction analysis cannot provide an unequivocal 
determination of the existence of r/-nucleus bound state. More direct mea- 
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surements are, therefore, necessary. 
PACS Numbers: 24.60.D, 13.75.G, 21.10.D 



lypeset using REVI^ 



I. INTRODUCTION 



The existence of //-mesic nucleus, a bound state of rj meson in a nucleus, was first 
predicted by us in 1986 The formation of such bound systems is caused by the attractive 
interaction between the r] meson and all the nucleons in the nucleus. The attractive nature 
of the interaction followed naturally from the work of Bhalerao and Liu [0 who found, from 
a detailed coupled-channel analysis of vrA^ ttN, ttN h-kN, and ttN rjN reactions, 
that near-threshold r]N interaction is attractive. We can easily understand this attraction 
by observing that the riN threshold is situated just below the iV*(1535) resonance. However, 
the attraction given in ref. [H is not strong enough to support a bound state in nuclei with 
mass number A < 12. This latter conclusion [||] was confirmed by Li et al. [§, who used 
a standard Green's function technique of many-body theory to study the formation of rj- 
mesic nuclei. Before proceeding further, we would like to mention that in the literature the 
?7-nucleus bound states are also called ?7-nucleus quasibound states. This is because these 
bound states have a finite width and they eventually decay. In this work, we shall simply 
call them bound states. 

If the existence of r^-mesic nucleus is experimentally confirmed, many new studies of nu- 
clear and particle physics will become possible. For example, because the binding energies 
of 7] meson depend strongly on the coupling between the rjN and the A^*(1535) channels 0, 
studies of ?7-mesic nucleus can yield additional information on the r]NN* coupling constant 
involving bound nucleons. Furthermore, the 77-mesic nuclear levels correspond to an ex- 
citation energy of ~ 540 MeV, to be compared with an excitation energy of ~ 200 MeV 
associated with the A- and E-hypernuclei. The existence of nuclear bound states with such 
high excitation energies provides the possibility of studying nuclear structure far from equi- 
librium. 

Several experiments 0,^0 using 7r+ beam, as motivated by the theoretical works of 
refs. J^flj were performed to search for r^-mesic nuclei. While these experiments could 
not confirm the existence of ?7-mesic nucleus, they did not rule out such possibility either. 
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More recently, Sokol et al. P,p!0| have claimed observing the ?7-mesic nucleus ^^C through 
measuring the invariant mass of correlated n^n pairs in a photo-mesonic reaction. Further 
confirmation of their experimental result, with improved statistics, is needed before arriving 
at a definite conclusion. 

The existence of these r^-mesic nuclear states depends on the value of the riN scatter- 
ing length a^AT. (For the scattering length, we use the sign convention of Goldberger and 



Watson [|TI1.) After 1990, many theoretical models were proposed for rjN interaction. Fits 



to various data have yielded very different a^^v i,|TljTl,[TlJTl,|T10,[T8||T|,|2|j2T|,|2|j2|]. The 



real part of a^Ar ranges from as low as 0.270 fm to as high as 1.000 fm while the 



imaginary part varies between 0.190 fm g] to 0.399 fm ||23|. This wide range of values of 
a^TV are summarized in table |. These very different values arise because a^^f is not directly 
measurable and must be inferred indirectly from other observables, such as the vrA^ phase 
shifts. In this latter case, the calculated a^^v depends strongly on the model used to relate 
the rjN to the ttN channels. As pointed out in ref. |T^, the inclusion of both the A^*(1535) 
and A^*(1650) resonances leads to larger scattering lengths. However, this general rule does 
not account for all the reasons. For example, in ref. ||T^ the inclusion of both these A^* 
resonances, but with fit to different reaction data, does not lead to a TZe arjN as high as 
^ 0.8 fm rather only to a TZe a.r^N ~ 0.6 fm. We caution the readers that it is premature 
to conclude which published o^at is the realistic one. Such a determination can be made 
only after r^-mesic nuclei are experimentally discovered, as only their experimental observa- 
tion can set a stringent constraint on the in-medium scattering length and enable one to 
differentiate between different theoretical models. 

Using final-state-interaction (FSI) analysis, several theoretical studies have suggested an 
O'-qN ~ (0.5 + 0.3?) fm and the possibility of the formation of bound ?7-nucleus states in light 



nuclei, such as ^'^He p^ , |25[] and deuteron p6| , |27|j28[| . We will revisit the FSI analysis in this 
work and show the limitation of the method. In spite of the fact that many predicted o^at 
(see table |) are much larger than the one used in the above-mentioned FSI analysis, to this 
date there has been no direct observation of these bound states in light nuclei. Consequently, 



we believe it is valuable to understand the current experimental situation by analyzing in 
detail the dynamics pertinent to the formation of ?7-mesic nucleus. 

We have learned from studying pion- nucleus scattering that, because of the Fermi motion 
and binding of the target nucleon, the pion-nucleon interaction in nuclei occurs at energies 
below its free-space value, i.e., below the energy that the pion-nucleon system would have 



when the nucleon is unbound As the above causes for the lowering of the interaction 
energy are also present in ?7-nucleus scattering, we believe that the ?7-nucleon interaction 
energy in r^-mesic nucleus formation will occur below its free-space threshold and this sub- 
threshold interaction can have important effects on the formation of //-nucleus bound states. 
This observation has motivated us to examine in this work the reliability of using the free- 
space TjN information, such as the rjN scattering length, to make predictions of the nuclear 
binding of the rj. 

We will, therefore, study the dependence of ?7-nucleus bound-state formation on off-shell 
behavior of rjN interaction in nuclei with A > 3 through comparing three different theoret- 
ical approaches. They are: (1) a fully off-shell microscopic optical potential that uses an 
off-shell rjN model; (2) a factorization approximation (FA) to the preceeding microscopic 
potential, where an energy shift parameter is used to calculate the rjN interaction at sub- 
threshold energies; and (3) an "on-shell" optical potential whose strength depends solely on 
the scattering length ar,N- The last approach has its root in the study of K- and vr-mesic 



atoms |]30| , pl| , p2| and is also extensively used in predicting 77-mesic nuclei. Finally, we men- 
tion that as the existence of 77-mesic nucleus has not yet been experimentally confirmed with 
certainty, there is no data to be fitted. Hence, all of our results as well as those of others 
are purely predictions. 

The paper is organized as follows. The three theoretical approaches are described in 
sec. ^ The important differences in their underlying reaction dynamics are outlined. Results 
and discussion are presented in sec. |T|. As we shall see, for a given rjN interaction model, the 
"on-shell" approach gives much stronger binding for 77 than does the microscopic calculation, 
making evident the important effects of off-shell dynamics. We will further show that these 



effects are of very general nature and are independent of tlie specific unitary off-shell model 
used in the comparative study. The ramifications of off-shell effects on both experimental and 
theoretical studies of bound states of rjA systems, particularly in systems with 3 < A < 12, 
are discussed. In view of the interest in the method of FSI, we also reexamine in more detail 
this method and its applications to the study of r7-^''^He systems. The conclusions of our 
study are summarized in sec. |V[ 



II. THEORY 

We calculate the binding energy and width of an ?7-nucleus bound system by solving 
the momentum-space relativistic three-dimensional integral equation 

^ V(k') + / dk<k' \V \ k> V'(k) = E^(k') , (1) 

using the inverse-iteration method Here < k' | | k > are momentum-space matrix 
elements of the ?7-nucleus optical potential V, with k and k' denoting, respectively, the initial 
and final //-nucleus relative momenta. The fi is the reduced mass of the //-nucleus system 
and E is the complex eigenenergy which we will denote as + = K^/2yU. For bound 

states, < and F^ < 0. As mentioned in sec. |, three different approaches to V are used 
by us to calculate E, and they are described below. 



A. Covariant //-nucleus optical potential 



In spite of its Schrodinger-like form, eq.(|I]) is covariant as it can be obtained from applying 



a specific covariant reduction to the relativistic bound-state equation \E' = G'oV\E'. The 
three-dimensional relativistic wave function and the covariant potential < k' | F | k > in 
eq.(|I|) are related to the fully relativistic ones by 



^(k) 



\i?(k 



(2) 



and 



< k' I y I k >= V^(k'2) < k' I V{W, k'\ r) I k > jR{k 



(3) 



In the above equations, 



(4) 



and 



M^ + Ma 



(5) 



E,(k) + E^(k) 

However, as a result of the apphcation of the covariant reduction, the zero-th components 
of the four-momenta k and k' are no longer independent variables but are constrained by 



= W - EA{k), k'^ = W-EAik') 



(6) 



The main advantage of working with a covariant theory is that the 77-nucleus interaction 
V can be related to the elementary rjN process by unambiguous kinematical transforma- 



tions |9 



The first-order microscopic //-nucleus optical potential has the form 

< k' I y I k > =J2jdQ<k', -(k' + Q) I t{^)^N^vN I k, -(k + Q) > 
x0;(-k'-Q)0,(-k-Q) , 



(7) 



where the off-shell rjN interaction t>riN^r]N is weighted by the product of the nuclear wave 
functions 4'*<f)j corresponding to having the nucleon j at the momenta — (k+Q) and — (k'+Q) 
before and after the collision, respectively. The is the rjN invariant mass and is equal 
to the total energy in the cm. frame of the rj and the nucleon j. It is given by |^ 



= [{W - Ec,{Q)}' - 



Mn + M, 



N- 



Q2 [M^ + Ma 



(8) 



where Q, Ecj and M^j are, respectively, the momentum, total energy, and mass of the 
core nucleus arising from removing a nucleon j of momentum —(k + Q) and binding energy 



I tj I from the target nucleus having the momentum — k. Equations (J^) and indicate 
that the calculation of V involves integration over the Fermi motion variable Q and requires 
knowledge of the basic t^^^riN at subthreshold energies. 

The matrix element of tr^N^riN in the //-nucleus system is related to the rjN scattering 
amplitude A in the rjN system by 

< k', -(k' + Q) I t{^),^N^nN I k, -(k + Q) > 



P > P) > (9) 



/E,(k')^7v(k' + Q)E^(k)E^(k + Q) 
where p and p' are the initial and final relative three-momenta in the riN cm. frame. We 
define the on-shell limit a.s p' = p = po and = Enipo) + Ej^{po) = ^/s^, where Po is the 
on-shell (asymptotic) momentum. A natural way of parameterizing A is 



^(V^, P', P) = ; -^(V^, P', P) , (10) 

An^^E,{p')EM{p')E,{p)Er,{p) 
so that {da/dQ)r]N^TiN =\ ^ The has the standard partial- wave expansion of a spin 
0-spin 1/2 system: 

^(V^,P',P) = tiT,2jAV^,P,P) + i^+^)tiT,2uiV^,p',p) ) Pii^) 



VPP I 

-ia-{px p') (4T,2i_ (v^, P'. P) - 4T,2i+ (v^. P'^ P) ) Peiz) ] , (11) 

where 2; = p ■ p', J± = £ ± 1/2, and T is the isospin of the rjN system and equals to 1/2. In 
the on-shell limit, 



^^(exp[2z52T,2,±(v^)]-lj • 



(12) 



fpfp 2ip 

The phase shifts 6^ are complex- valued because the thresholds for t]N ttN and t]N ttttN 
reactions are lower than the threshold for r]N scattering. When p ^ 0, 5^ p^^+^a^^) and 



Vpp 



P''a^l2,, . (13) 



The a2T^2j and 037! 2^ are, respectively, the (complex) riN scattering length and volume. Near 
the threshold, only the s-wave term, t^i in eq.(|TT]), is important. 



Different off-shell models give different off-shell extensions of A to kinematic regions 
where p p' and y/s] 7^ ^/s^ . In the separable model of ref. 0, the off-shell amplitude is 
given by 



ta{y/s],p,p) = K{^,p,p) ^p'p \^ D {^) 



(14) 



with 



K= -^JE^{p')Er,{p')E,{p)EM , (15) 



Na = hai^/s],p)hai^,p) OC |^ {p' pY V t{p')v ^{p) , (16) 

and 

= ^ - - s;;(y5-) - s°(y5-) - KAV^) ■ (i7) 

Here a is a short-hand notation for the quantum numbers (£, 2T, 2j) of the isobar resonance 
a. The is the bare mass of the isobar a and S^, S", and S"^ in eq. ([T7D are the self- 
energies of the isobar a associated, respectively, with its coupling to the rjN, ttN, and tcttN 
channels [@. The coupling constants and form factors are denoted by g and v. At the 
rjN threshold, only the s-wave 77 interaction is important, which limits the isobar a to 
A^*(1535) or a = {i,2T,2j) = (0, 1, 1). Clearly, different models will have different off-shell 



extensions in energy and momenta for ta{y/Sj,p' ,p). However, they should all satisfy eq.(|T2D 
in the on-shell limit. 

B. Factorization approximation 

We define the factorization approximation (FA) by taking the rjN scattering amplitude 
in eq.(^ out of the Q integration at an ad-hoc fixed momentum < Q >: 

< k' I VpA I k > = < k', -(k'+ < Q >) I t(v^)^7v-.^jv I k, -(k+ < Q >) > 

X /(k' - k) , (18) 



where 

/(k' ~k)=^JdQ 0*(-k' - Q)0,(-k - Q) , (19) 

j 

is the nuclear form factor having the normahzation /(O) = A. In eg . (p!8D , tr^N^rjN is still 
defined by the same functional dependences on various momenta and energies as given by 
eq.(|l^, except that p' and p are now determined from an ad-hoc momentum < Q > in the 
77-nucleus system and that the interaction is given by an ad-hoc energy i/f • The choice of 
< Q > is not unique. One option is to take an average of two geometries corresponding, 
respectively, to having a motionless target nucleon fixed before and after the t]N interaction. 
This leads to 

< Q >= - (^) (k' - k) . (20) 

This choice has the virtue of preserving the symmetry of the t-matrix with respect to the 
interchange of k and k'. (There are other possible schemes; see for example ref. @].) An 
inspection of eq.(||) suggests that it is reasonable to assume a/I = Mr^ + — A = ^/s^— A, 
with A being an energy shift parameter. In tcN scattering, the downward shift A that fit 
the data was determined to be ~ 30 MeV 



C. "On-shell" optical potential 

We will use the term "on-shell" optical potential for the optical potential where on- 
shell hadron-nucleon {hN) scattering length is used to generate the hN interaction. In 
the literature, the first-order low-energy hadron-nucleus "on-shell" optical potential is often 



given as ||36| , |37 



X 



y I k' I I k I a) ^ , . 



where is the hadron mass, /i the hadron-nucleus reduced mass, / the nuclear form factor 

''r]N 

10 



normalized as /(O) = A. For ?7-mesic nuclei, Mh = and a^f'^f is the s-wave t]N scattering 



length. As has been shown in ref. U corresponds to Vfa with no energy shift (A = 0) 
and in a static hmit of the target nucleon. In this latter limit, the r]N relative momenta are 
p = k/(l + M^/Mn) and p' = k!/{l + M^/Mn). In addition, k' • k = p' • p. 

III. RESULTS AND DISCUSSION 

The binding energies and half-widths of r^-mesic nuclei given by the off-shell microscopic 
calculation are presented in table |T|. The covariant optical potential [eq.(^] used in the 
calculation is based on the model of ref. 0. We use this model to demonstrate the effect of 
subthreshold rjN interaction on the formation of r^-mesic nucleus. The solutions are obtained 
with the rjN interaction parameters QriNa, Ma, and A^tvq determined from the vrA^ phase 
shifts of Arndt et al. The p-wave and d-wave interactions are also attractive at the threshold 
but their magnitudes are very small and have negligible effect on and F^. The nuclear 
wave functions in eq.(|^) are derived from the experimental form factors with the proton 



finite size corrected for ||3^. Details of the calculation can be found in refs. As can 

be seen from the table, the binding energy increases as the nucleus becomes heavier. In 
addition, the number of nuclear orbital in which the t] is bound (Is, 2s, Ip, etc.) increases 
with increasing mass number A. The reason for this trend is discussed in ref. We would 
like to point out that our microscopic calculation does not use on-shell scattering length as 
an input. 

We emphasize that the reason we use the off-shell model of ref. in our comparative 
analysis is because we have at our disposal the detailed unitary off-shell momenta and energy 
dependences of that model, which allow us to study the off-shell effects in the formation of 
77-nucleus bound-state. Such off-shell information is not readily obtainable from other r]N 
models, either because the models cannot be extended to off-shell domain or because of the 
elaborate computation required to generate the needed off-shell p,p', and -^i dependences 
of the riN interaction. Fortunately, as we shall see later, the effects of off-shell dynamics on 
any reasonable off-shell model will be qualitatively similar. 
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In table |T|, we present the bound-state solutions obtained from using the factorized 
covariant potential Vfa [eq.(p!8D] with A = 0, 10, 20, 30 MeV. The same interaction 
parameters were used. The nuclear form factors p9| , |i0[] used in the calculations are sum- 
marized in table A comparison between tables ^ and ^ indicates that the FA results 
with A = 30 MeV are very close to the off-shell results. This value of A is similar to the 
one found in pion-nucleus elastic scattering p3| and can be understood by noting that the 



average nuclear binding and Fermi motion amount to about 30 MeV downward shift p9| of 
the hadron-nucleon interaction energy ^/s. Our full off-shell dynamical calculations indicate, 
therefore, that the rjN interaction in rj bound state formation takes place at energies about 
30 MeV below the (free-space) threshold. 

The subthreshold nature of the hadron-nucleon interaction in a nucleus is also evident in 
i^-mesic and vr-mesic atoms. For example, although the free-space KN scattering length is 
repulsive, the effective KN scattering length needed to fit the i^-mesic atom data is attrac- 
tive [^]. The sign change of the scattering length can be easily understood. Indeed, one 
may note that the A (1405) resonance is situated about 26 MeV below the KN threshold 
||30|| . Using a downward shift of 30 MeV found in this work, we can see that the in- medium 
KN interaction occurs actually below the resonance. Hence, an attraction arises (see sec. |). 
In the TT-mesic atom studies, a strong s-wave repulsion of the effective vrA^ interaction is in- 
dicated by the data. This is in sharp contrast to the free-space vrA^ s-wave interaction where 
there is a near-perfect cancellation between the 5*31 and Su scattering lengths. Bhalerao 
and Shakin showed that the strong s-wave repulsion is due to different energy depen- 
dences of the S'31 and Sn interactions, which make the cancellation no longer near complete 
at subthreshold energies. The mesic-atom experiments, therefore, indicate unambiguously 
that in bound-state formation the hadron-nucleon interaction inside a nucleus occurs at 
subthreshold energies. 

The results of the "on-shell" optical potential [eg. (pi])] are given in table |V] for two 
different values of a^Tv- The first one is the same one used for tables |I| and |T|, the other 
is chosen from table |, such that it has an imaginary part very similar to the former one. 
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The nuclear form factors are the same as those used in the FA calculations. For these two 
scattering lengths, no bound state can exist in ^He. Upon comparing the third column of 
table [V| with the off-shell calculation (table |I|), we can see that the "on-shell" approximation 
predicts more strongly bound //-mesic nuclei. Also, as expected, the "on-shell" results for 
O'rjN = (0.28 + 0.19i) fm are similar to those of FA with A = MeV. We further note from 
table |V| that the effect on decreasing the calculated nuclear binding energies as caused by 
the decrease in TZe ajjN is not a linear function of the nuclear mass number A. For light 
nuclei (A < 12) a decrease of TZe a^jN by a factor of ~ 2 causes to decrease by a factor 
greater than 10. 

While the in-medium rjN interaction strength is a function of the off-shell variables 
p'jP^ \fsi can appreciate the main feature of the off-shell effects by considering only the 
a/s, or the energy, dependence. Our study of the three different optical potentials indicates 
that this strength decreases with energy. However, specific subthreshold energy dependence 
is model-dependent. We may define a phenomenological reduction factor R by 

^N^rjNiy/s^- ^,P',P) = RiA)trjN^r,Ni^/s^,P',P) , (22) 

and introduce an effective in-medium scattering length by 

a;^J = i?(A) . (23) 

In the model ofref. §, i?(A) = ^ D(^)/[(y5^- A)D(^ - A)] [eqs.(|li and (|ll- 
(P)]. Our calculation indicates that a/J = (0.26 + 0.13i), (0.25 + O.lli), (0.23 + 0.09?) fm 
for A = 10, 20, 30 MeV, respectively, while at the threshold a^jv = (0.28 + O.lQi) fm. The 
R{A) is about 0.89 at A = 20 MeV and 0.82 at A = 30 MeV. This reduction is the origin of 
the A-dependence of the and in table |ITT[ For the scattering length of the model of ref. 
||12|| , it was mentioned that at 20 and 30 MeV below the threshold, a^^J = (0.49 -|- O.lOz) fm 
and (0. 45-1-0. 08i) fm, respectively. Upon comparing these values with a^^v = (0. 75-1-0. 27i) fm 
given by that model at the threshold, we see a reduction of more than 1/3 of the real part 
of arjN (i-e., R = 0.6). It seems that higher is the TZe ar,N at the threshold greater is the 
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subthreshold reduction. We note that at 20 to 30 MeV below the threshold, both models 
give an imaginary part of a^Jj^ about 0.09i fm. This is because the Xm a^-^ is related to 
the reaction channel. At the subthreshold region, the only reaction channels that the rjN 
channel can couple to are the vrA^ and vrvrA^ channels which were taken into account by 
both the models P,p!^. On the other hand, the real part of the effective scattering length 
is still very model dependent, though in lesser degree than at the threshold. Because of 
this model dependence, it is not possible to guess the subthreshold reduction for the other 
models listed in table | for which the off-shell dependence cannot be easily reconstructed 
from the corresponding publications. We conclude from the two models analyzed above that 
a substantial reduction of attractive strength, IZe arjN, must occur at subthreshold energies. 

This reduction in attractive strength is of very general nature as it is a direct consequence 
of the A^*(1535) resonance. We recall that the attraction between the t] and nucleon at low 
energies arises because the r]N threshold is situated below this resonance. However, as the 
energy shift A becomes larger, the rjN interaction energy moves farther downward away 
from the resonance. Any reasonable off-shell model leads, therefore, to a reduced attrac- 
tion. Consequently, calculations making use of free-space scattering length (corresponding 
to A = 0) necessarily overestimate e^. Because of the large sensitivity of the binding ener- 
gies to ?7- nucleon interaction in light nuclei and because Tie a^{J <^ TZe a^^r, many //-nucleus 
bound states in very light nuclei, as predicted by using some of the a^^v in table |, may not 
exist in real situation. We caution, therefore, against using "on-shell" (or scattering-length) 
prediction as guide for searching ?7-mesic nuclei. 

In view of the existing interest in searching for ?7-nucleus bound states in light nuclear 
systems, we believe that it is informative to determine values of the "minimal" scattering 
length, a^^, which represents the least value of an a'^-^J that can bind the r] into an Is nuclear 
orbital. Clearly, the real and imaginary parts of this scattering length are not independent of 
each other. We fixed Xm a^" to 0.09 fm, as suggested by the two off-shell models discussed 
above, and searched for TZe aj^". The results for several light nuclei are given in table [VT| . 

It is interesting to revisit the FSI analysis and to examine the results of refs. P^J^ 

14 



in the light of our findings. The FSI analysis was first applied by Wilkin ||2^ to rj-^He 
system which is a final state of the pd rj'^B.e reaction. The analysis raised a great deal 
of theoretical interest and was later extended to the study of rj-^He system |^^, which 
is a final state of the dd —>■ rj'^B.e reaction. Let us denote the 77-nucleus scattering length as 



Aq = Af> + iAi. According to the Watson FSI theory for a weak transition and a strong 
FSI, one can approximate the total reaction amplitude, /, at low energies by 



/ 



(24) 



I |2 {l + Ajpr,)^ + {Anp^y ' 

where Pr^ is the cm. momentum of the t] meson. The fs is the transition amplitude and 



was treated as a constant in ref. ||2J]. The unitarity condition requires Aj > 0. Hence, 
I / P increases monotonically when decreases. A good fit to the p^-dependence of | / p 
was obtained by Wilkin who used an a^iv = (0.55 ± 0.20 + O.SOi) fm, which led to an 
y4o(?7^He) = (—2.31 + 2.57i) fm |2^. In a later work, Wilkin and his collaborators |^ 



used a^AT ^ (0.52 + 0.25i) fm which gave Ao{r]^Re)^ (—2.3 + 3.2i) fm and Aolrj'^Re) a; 
(—2.2 + l.lz) fm, respectively. Upon introducing these Aq into eg. (p^ , they obtained a good 
representation of the experimental p^j-dependence of | / p for both '^'^He. Because in the 
case of a real potential, a negative scattering length has the possibility of corresponding 
either to a repulsive interaction or to an attractive interaction that supports a s-wave bound 
state the findings of ref. |^ has raised the hope that ?7-nucleus bound states might 
exist in ^'^He. 

A closer look at the value of the r^-helium scattering lengths can provide partial answer 
to the question as to whether rj can be bound in ^''^He nuclei. This is because for a complex 
potential the existence of a bound state imposes a constraint on the complex scattering 
length. Equation ( p4D indicates that, if a weakly-bound state exists, the amplitude / will 
have a pole, Ppoi, in the complex p^-plane at 

—i —Aj — iAr 
Vvoi = ^ = ^ , (25) 

where d = A\ + A?^ and Aj is positive. The condition for a bound state is IZe p^^i < 0. This 
requires that | Ai \<\ Ar |. It is easy to see that the AoijfYie) of ref. does not satisfy 
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this inequality while the Ao{r]'^B.e) does. Our optical-potential calculations have verified 



this, namely, with the o^at used in ref. |^2|] we have found that there is no bound state in 
^He but there is one in ^He. In other words, A/j < is not a sufficient condition for having 
a bound state. Consequently, in spite of the much greater slope of | / p with respect to pr. 



in the 77- '^He system , we conclude that there is no ?7-nuclear bound state in ^He. For the 
same reason, 77 cannot be bound onto the deuteron because none of the various calculated 
AQ{r]d) of ref. can satisfy the condition | Aj \<\ Aji \. On the other hand, the formation 



of a bound state in ^He remains a possibility. 

This possibility is, however, hampered by the fact that | / p is insensitive to the sign 
of Aji, as can be seen from eq. . Indeed, we have found an Aq that has a positive real 
part and can equally describe the data. For example, with a^^r = (0.30 + 0.09i) fm, we have 
obtained from our optical potential, eq.(pT[), an Ao(?7^He)= 2.10 + 2.88i fm which can well 
describe the '^He data. Similarly, a good representation of the ^He data can be obtained 
with an y4o(?7^He) = (0.80 + 1.75z) fm by using a^^v = (0.16 + 0.09z) fm. The quality of 
representation is shown in fig.|l]. In both cases A^ > which corresponds to the situation 
where the interaction is attractive but not strong enough to support a bound state. Hence, 
we have shown that FSI can also admit solutions corresponding to having no 77 bound state. 
We emphasize that our above analysis does not imply our preference for weaker ar^N but 
rather to exemplify numerically that FSI cannot provide a unique answer concerning the 
formation of bound state. 

At this point, several comments are in order. Firstly, although we need different a^^v for 
^He and '^He while in ref. only one arjN was needed, our result is reasonable because it 



reffects the fact that the transition amplitudes fs in the two reactions may not be treated 
as p^-independent and, thus, do not scale each other by a multiplicative constant. Indeed, 
by using the two-step model of ref. |^^, Willis et al. ||^ cannot fit the data of | f{P'q) P for 



the ?7-^He and 77-^He systems with the same Qr^N- In addition, one may note that second- 
order ?7-nucleus optical potential is important in light nuclei. This potential is sensitive to 
the nucleon-nucleon (NN) correlation which is very different in ^He and "^He, as evidenced 
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by the fact that the matter densities derived from the measured charge densities of these 
two nuclei have spatial dependences that do not scale each other. Inclusion of second-order 
optical potential will lead to different modifications of y4o(?7^'^He). Consequently, when only 
first-order optical potential is used to generate the Aq in fitting the data, the effective o^at 
will be different. In a similar manner, the different A^A^ correlations in the two helium nuclei 
also contribute to the nonscaling of the fs- Indeed, many reaction processes contribute to 
the transition matrix element. An important subset of the processes are those that cannot be 
approximated by the pp ir^d and pd tt'^ ^He doorway mechanisms @]. This is because 



in these processes the intermediate pion interacts with two different nucleons that do not 
belong to the same deuteron cluster in the nucleus. Hence, these processes depend strongly 
on A^A^ correlations having momentum dependence that cannot be obtained from using the 
deuteron wave function. The non-simple nature of Jb is further exemplified by the work of 
Santra and Jain |^5[, who can describe the ?7-'^He data even without FSI when the exchanges 
of p and other mesons are included in the calculation of Jb- While we strongly believe that 
FSI must be taken into account, ref. |^5| does indicate that there are many more aspects 



of the dynamics left to be thoroughly investigated before a definitive conclusion could be 
drawn from the ?7-^'^He data. 

Secondly, our use of Xm a'^^J = 0.09 fm is merely suggested by the 30-MeV downward shift 
of the effective r]N interaction energy, as has been discussed after eq.(^). With Xm a'^^J= 
0.09 fm, the helium data constrain the real part of the effective scattering length to be 
between 0.16 and 0.30 fm. However, we emphasize that the value of 0.09 fm is model- 
dependent and it is important to see if different riN off-shell models would give a similar 
subthreshold value. Consequently, our finding of repulsive scattering length for the helium 
systems should not be interpreted as if Xm a^;j^=0.09 fm were the final answer, but rather 
as an example showing the nonuniqueness of the solution to the problem. In view of the 
many needed improvements in the modeling of the transition amplitude, searching a same 
repulsive effective o^at for both data sets is beyond the scope of this work. 

One would certainly like to be able to use the ?7-^'^He data to obtain a more stringent 
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constraint on the effective subthreshold rjN interaction. However, we conclude from the 
above discussion that the currently used FSI analysis cannot provide an unequivocal answer 
not only because there is no simple correspondence between the sign of and the existence 
of an r^-nucleus bound state, but also because the assumption of a p^-independent fs is 
uncertain. Consequently, our main message is that one should not rely on using on-shell 
scattering length and that direct detection of bound states is necessary. 

Finally, we have also examined effects of nuclear form factor on the binding of rj. The 
results for ^He are presented in table [VI 1| . An inspection of this last table indicates that 
the binding of t] is not very sensitive to the finer aspect of a realistic nuclear form factor. 
A similar situation has also been noted for ^He. The binding energies and widths of the 
rj as given by using the 3-parameter Fermi or the Frosch form factor are close to each 
other. The most important effect on the formation of r^-mesic nucleus is, therefore, from the 
subthreshold dynamics of the ?7N interaction. 



IV. CONCLUSIONS 

Our study shows that calculated binding energies and widths of 77-nucleus bound states 
strongly depend on the subthreshold dynamics of the t]N interaction. The results of mesic 
atoms and the present analysis indicate that the average rjN interaction energy in mesic- 
nucleus formation is below the threshold. What matters for the bound-state formation is 
not the 7]N interaction at the threshold but the effective in-medium interaction. Because 
the subthreshold behavior of riN interaction is very model dependent, we believe that it 
is useful for theorists to publish not only the ?7-nucleon scattering length a^Tv, but also 
the corresponding subthreshold values as a function of A. Before the availability of this 
information, we suggest to look for bound states in nuclear systems much heavier than 
those indicated by on-shell scattering length. 

The FSI analysis represents an interesting approach. Since the analysis by itself cannot 
provide a definitive answer as to whether there is an attractive ?7-nucleus interaction strong 
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enough to bind the rj, more direct measurements such as the one used in ref. |^ are necessary. 

The downward shift in the effective interaction energy can lead to a substantial reduc- 
tion of the attraction of in-medium r^-nucleon interaction with respect to its free-space value. 
Consequently, predictions based upon using free-space rjN scattering length inevitably over- 
estimates the binding of 77. This overestimation of the binding, as revealed by this study, 
has never been taken into account in discussing r^-nuclear bound states. One must bear 
this in mind when using the predictions given by such calculations as guide in searching for 
77-nucleus bound states. 
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TABLES 

TABLE L Eta- nucleoli s-wave scattering lengths a^N- 



a^AT (fm) 



Formalism /Reaction 



Reference 



0.270 + 0.220i 
0.280 + 0.190i 
0.281 + 0.360i 
0.430 + 0.394i 
0.579 + 0.399i 
0.476 + 0.279i 
0.500 + 0.330i 
0.510 + 0.210i 
0.550 + 0.300i 
0.620 + 0.300i 
0.680 + 0.240i 
0.750 + 0.270i 
0.870 + 0.270i 
1.050 + 0.270i 
0.404 + 0.343i 
0.876 + 0.274i 
0.886 + 0.274i 
0.968 + 0.281i 
0.980 + 0.370i 



Isobar model 
Isobar model 
Photoproduction of r] 



Electroproduction of rj 
pd — >^He er] 
Isobar model 

Coupled T-matrices 
Effective Lagrangian 
Coupled A'-matrices 
Coupled A-matrices 

Coupled T-matrices 



Bhalerao and Liu 

ibid 

Krusche 
ibid 
ibid 

Tiator et al. |2S 



Wilkin |24| 



Sauermann et al. |14] 
ibid 

Abaev and Nefkens 116 
Kaiser et al. |jl 
Green and Wycech [12 
Green and Wycech [13 
ibid 



Batinic et al. [18] 
Batinic and Svarc [19| 
ibid 

Batinic et al. 



Coupled T-matrices 



Arima et al. [21 



23 



TABLE II. Binding energies and half-widths (both in MeV) of ry-mesic nuclei given by the full 
off-shell calculation. The solutions were obtained with the rjN interaction parameters determined 
from the ttN phase shifts of Arndt et al. No bound state solutions of eq.(l) were found for ^4 < 12. 



Nucleus 


Orbital {ni) 




12c 


Is 


-(1.19-F3.67i) 


160 


Is 


-(3.45 + 5.38z) 


26Mg 


Is 


-(6.39 + 6.60i) 


40Ca 


Is 


-(8.91 -F6.80i) 


90Zr 


Is 


-(14.80 + 8.87z) 




Ip 


-(4.75 + 6.70z) 


208p|j 


Is 


-(18.46 + lO.lH) 




2s 


-{2.37 + 5.82i) 




Ip 


-(12.28 + 9.28z) 




Id 


-(3.99-F6.90i) 
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TABLE III. Binding energies and half-widths (both in MeV) of 77-mesic nuclei obtained with 
the factorization approach for different values of the energy shift parameter A (in MeV) . 



Nucleus Orbital {n£) 


A = 


A = 10 


A = 20 


A = 30 


12c 


Is 


-(2.18 + 9.96z) 


-(1.80 + 6.80i) 


-(1.42 + 5.19z) 


-(1.10 + 4.10z) 


160 


Is 


-(4.61 + 11.57?;) 


-(3.92 + 8.13i) 


-(3.33 + 6.37i) 


-(2.84 + 5.17z) 


26Mg 


Is 


-(10.21 + 15.41i) 


-(8.95 + 11.17z) 


-(7.94 + 8.97z) 


-(7.11 + 7.46z) 


40Ca 


Is 


-(14.34 + 17.06z) 


-(12.75 + 12.55i) 


-(11.53 + 10.2H) 


-(10.51 + 8.59z) 


902r 


Is 


-(21.32 + 18.59i) 


-(19.15 + 13.97i) 


-(17.58 + 11.54i) 


-(16.29 + 9.84i) 




Ip 


-(8.27 + 16.0H) 


-(7.19 + 11.47i) 


-(6.23 + 9.48i) 


-(5.40 + 7.94i) 


208 p]^ 


Is 


-(24.06 + 19.18z) 


-(21.88 + 14.44z) 


-(20.28 + 11.96) 


-(18.96 + 10.22z) 




2s 


-(4.89 + 11.04i) 


-(3.67 + 8.28i) 


-(2.81 + 6.79i) 


-(2.12 + 5.72i) 




Ip 


-(18.33 + 18.97z) 


-(16.31 + 14.27z) 


-(14.81 + 11.79z) 


-(13.56 + 10.060 




Id 


-(8.27 + 14.07z) 


-(6.17 + 10.56z) 


-(5.58 + 8.7H) 


-(4.66 + 7.4H) 
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TABLE IV. Nuclear form factors used in the factorization approach and scattering-length 
approach. 



Nucleus 



Form factor 



Parameters^ 



•^He 

^He 

6Li 
9Be 

12c 
160 
26Mg 

40Ca 
902r 

208 p]^ 



Hollow Exponential 
Gaussian 

3-parameter Fermi 
Frosch model 
Modified Harmonic Well 
Harmonic Well 
Harmonic Well 
Harmonic Well 
Harmonic Well 
Harmonic Well 
2-parameter Fermi 

2- parameter Fermi 

3- parameter Gaussian 
2-parameter Fermi 



a = 1.82 fm 
a = 1.77 fm 

c = 1.01 fm, z = 0.327 fm, w = 0.445 fm 
a = 0.316 fm, b = 0.680 fm 
ai = 1.71 fm, 02 = 2.08 fm 
a = 2/3, a = 2.42 fm 
1, a = 2.45 fm 
1, a = 2.42 fm 
a = 4/3, a = 2.53 fm 
a = 1.6, a = 2.75 fm 
c = 3.050 fm, z = 0.524 fm 
c = 3.510 fm, z = 0.563 fm 
c = 4.500 fm, z = 2.530 fm, w = 0.20 fm 
c = 6.624 fm, z = 0.549 fm 



a 



a 



Ref. for A = 3 — IQ and ref. |^0| for the rest of the nuclei 
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TABLE V. Binding energies and half-widths (both in MeV) of 77-mesic nuclei given by the 
scattering-length approach for two different values of the scattering length a^iv. A blank entry 
indicates the absence of bound state. No bound state exists in ^He. 



Nucleus Orbital {nt) o^at = (0.28 + 0.19z) fm o^at = (0.51 + 0.2li) fm 



^He^ 


Is 




-(6.30 + 11.47i) 


^Li 


Is 




-(3.47-|-6.79z) 


9Be 


Is 




-(13.78 -1- 12.45i) 




Is 


-(0.93 + 8.70) 


-(15.85 + 13.05i) 


iiB 


Is 


-(2.71 + 10.9H) 


-(20.78 -M5.42?) 


12c 


Is 


-(2.91 + 10.22z) 


-(19.61 + 14.20?) 


160 


Is 


-(5.42 + 11.43i) 


-(23.26 + 14.86i) 




Ip 

IT 




-(0.95-F7.72z) 


26Mg 


Is 


-(11.24-1- 14.76z) 


-(33.11 -1- 17.73i) 




\p 

Ir 




-(13.41 + 12.33i) 


40Ca 


Is 


-(15.46 -1- 16.66i) 


-(38.85 -1- 19.16i) 




2s 




-(5.59-F6.14z) 




\p 


— (1.22 -1- 10.58^) 


— (22.84 + 14.32i) 




Id 




-(4.28-K9.52i) 


902r 


Is 


-(22.41 + 19.97i) 


-(48.40 + 22.60z) 




2s 




-(26.07 + 10.07i) 




\p 


-(10.18 -M4.33i) 


-(31.53 -M5.93i) 




2p 




-(18.51 + 8.57z) 


208pb 


Is 


-(24.55 + 19.57i) 


-(50.27 + 21. 42i) 




2s 


-(10.56 -M3.32i) 


-(22.27 + 11.50i) 




Ip 


-(20.19 + 19.05z) 


-(34.03 + 10.03i) 




2p 




-(1.89 + 3.75i) 
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Id -(12.22 + 16.07i) -(27.89 + 12.17i) 



^ Form factor used is the 3-parameter Fermi. 



TABLE VI. Values of a^^" for nuclei with mass number A < 10. The Im a"^}^ is fixed at 0.09 
fm (see the text). 



Nucleus Nuclear form factor a^ff' (fm) 

^He Hollow exponential 0.49 + 0.09i 

^He 3-parameter Fermi 0.35 + 0.09i 

•^Li Modified harmonic weU 0.35 + 0.09z 

^Be Harmonic Well 0.24 + 0.09z 



TABLE VH. Dependence of binding energies and half-widths (both in McV) of ry-'^Hc bound 
state on two different form factors (hollow exponential and Gaussian) for a few values of the 
scattering length a^jv- All bound states are in the Is orbital. 



a-qN fm Hollow Exponential Gaussian 

0.680 0.240i -(3.74 7.89i) -(4.16 7.93i) 

0.750 + 0.270i -(6.24 + 9.94i) -(6.79 + 9.94z) 

0.876 0.274i -(11.86 11.26z) -(12.59 -Fll.lOi) 
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FIGURES 




0.1 0.2 0.3 

[1/fm] 

FIG. 1. Dependences of | / p on p^. Data for the ?7^He (crosses) and r/^He (open circles) 
systems are from refs. and |42|, respectively. Solid curves are the results of using the scattering 
lengths having Ar > 0. 
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